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Abstract 

In the seesaw model with hierarchical Dirac masses, the neutrino mixing an- 
gle exhibits the behavior of a narrow resonance. In general, the angle is strongly 
suppressed, but it can be maximal for special parameter values. We delineate the 
small regions in which this happens, for the two flavor problem. On the other hand, 
the physical neutrino masses are hierarchical, in general, except in a large part of 
the region in which the mixing angle is sizable, where they are nearly degenerate. 
Our general analysis is also applicable to the RGE of neutrino mass matrix, where 
we find analytic solutions for the running of physical parameters, in addition to a 
complex RGE invariant relating them. It is also shown that, if one mixing angle is 
small, the three neutrino problem reduces to two, two flavor problems. 
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1 Introduction 



The exciting development of recent experiments |T[ has offered strong evidence for the 
existence of neutrino oscillation, from which one can infer about the intrinsic properties 
of the neutrinos. While the neutrino masses (mass differences) are found to be very tiny, 
there is a major surprise for the mixing angles. It is found that at least one, and possibly 
two, of the three mixing angles are large, or even maximal. This is in stark contrast to 
the situation in the quark sector, where all mixing angles are small. 

Theoretically, the seesaw model [@ is very appealing in that it can offer a natural 
mechanism which yields small neutrino masses. However, owing to its complex matrix 
structure, it is not obvious what the implied patterns of neutrino mixing are. In a previous 
paper we found a parametrization which enabled us to obtain an exact solution to 
the two flavor seesaw model. When one makes the usual assumption that the Dirac mass 
matrix has a strong hierarchy, the physical neutrino mixing angle exhibits the narrow 
resonance behavior. For generic parameters in the Majorana mass matrix, the physical 
neutrino mixing angle is strongly suppressed. However, if the parameters happen to lie 
in a very narrow region, the mixing can be maximal. 

In this paper we will expand on our earlier investigations and discuss in detail the 
behavior of the neutrino mixing matrix in the seesaw model. As was shown before, if we 
assume that the Dirac mass hierarchy is similar to that of the quarks, the problem has 
three relevant parameters associated with the Majorana sector, namely, the mixing angle, 
the ratio of masses, and their relative phase. We will present plots of the physical neutrino 
mixing angle and their mass ratio in the 3D parameter space. These will offer a bird's-eye 
view of their behaviors. In particular, the neutrino mixing angle is only appreciable in a 
very small region, which we exhibit explicitly. Furthermore, this region is complementary 
to the region in which there is appreciable physical neutrino mass hierarchy. Thus, roughly 
speaking, the seesaw model divides the 3D parameter space in two parts. There is a very 
small region in which the mixing angle is large, at the same time the neutrino masses are 
nearly degenerate. For most parameters, the mixing angle is small but there is a strong 
hierarchy in the mass eigenvalues. An exception to this picture is when the Majorana 
matrix has extreme hierarchy and very small mixing angle. In this tiny region, the physical 
neutrinos can be hierarchical and simultaneously their mixing angle is large. 

The solution to the seesaw problem is most transparent in the parametrization in- 
troduced in Ref. 0. However, it is useful to make connections with the usual notation, 
where individual matrix elements are regarded as independent parameters. We obtain 
relations which clarify the roles played by the various parameters. They enable one to 
gain insights in understanding the numerical results presented in the 3D plots. 

The general analysis of symmetric and complex matrices turns out to be useful in other 
applications. Our method can be used to yield an analytic solution of the renormalization 
group equation (RGE) of the effective neutrino mass matrix. In addition, we obtain a 
(complex) RGE invariant which relates the running of the mixing angle and the complex 
mass ratio. The detailed analysis of the structure of the seesaw matrix also suggests a 
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universal picture for the quark as well as the neutrino mass matrices. While the quarks 
have generally small mixing angles and hierarchical mass ratios, if one assumes that the 
Majorana matrix itself is of the seesaw type, the effective neutrino mixing angle can be 
naturally large. 

Finally, we turn to a discussion of the three neutrino problem. Although the principle 
involved here is the same as in the two neutrino problem, the algebra with the Gell-Mann 
A matrices is far more complicated than that of the Pauli a matrices. We are unable 
to obtain a general solution in this case. However, it is quite well-established that one 
of the neutrino mixing angles is small IQ. In this case, an approximate solution can be 
obtained. It turns out that, to lowest order, the three neutrino problem can be reduced 
to two, two-flavor problems. This solution can thus accommodate the "single-maximal" 
or "bimaximal" solutions that have been considered in the literature. 



2 The two flavor problem 

In a previous paper 0, 0, an exact solution was obtained for the two flavor seesaw model. 
In this section, in addition to a summary of the earlier paper, further results will be 
presented. 

For two flavors, the seesaw model 



rrii, = m^Mj^^nr^ 



can be written in the form 
m.. = u( 



(1) 



(2) 



Let us introduce the parametrization 
I nil \ 



m2 



^Jm^rn^e ^"^ , ^ = - ln(m2/mi); 



(3) 



Ml 



-1 



M1M2 



■ln(M2/Mi); 



V, 



R 



-«/3o"2gi7o"3 



(4) 
(5) 

Thus, in the basis in which rriD is diagonal, P is the mixing angle for M^^ while ±27 
are the phases of the eigenvalues. This parametrization shows clearly that the relevant 
variables in the diagonalization of m^, are ^, a, P, 7 and 77. Of these, it is usually assumed 
that 1712/ mi can be identifled with the known quark mass ratio. Also, a can be absorbed 
into U as part of the phase of the Dirac mass eigenvalues. For [/ ~ /, in particular, it 
becomes the phase of the charged leptons and is not observable. 
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Note also that, apart from an overall constant, m,^ is a product of 2 x 2, complex 
matrices with det = +1, i.e., it is an element of SL{2, C). Thus, we can identify m,^ with 
an element of the Lorentz group, with and r] interpreted as rapidity variables. 

To find the effective neutrino mixing matrix, we need to rearrange the matrices in 
in a different order 

^ gia.'a3g-iea2gi0<T3^ (7) 

u' = uj + a, A = ^ ln(/i2/Ati)- (8) 

Here, the physical neutrino masses are given by /ii and /i2, with their ratio given in terms 
of A, while 40 is their relative phase. We have also absorbed the phase a into uj' . The 
physical neutrino mixing matrix is given by UW , so that W is the induced mixing matrix 
from the seesaw mechanism. The left-handed Dirac mixing, f/, in analogy to the quark 
sector, is often taken to be close to the identity, f/ ^ /. In the following we will concentrate 
on the behavior of W only, corresponding to U ~ /. However, when necessary, U can 
always be included in the final result. 

As was shown before, the solution for W corresponds to that of the velocity addition 
problem in relativity, and one can readily obtain the answer by manipulating the Pauli 
matrices. We have 

tan2cj = ; — -, (9) 

Sjj coth2^-cos2/3' ^ ' 

tan2^ = sin2/3/(cos2^cosh20 ^^^^ 

cos 2(3 — Tjji tanh 2^ — S/ tan 2uj 

cosh 2 A = cosh 2^ cosh 2fj — cos 2/3 sinh 2^ sinh 2f/, (11) 
where X = \ + icp, ^ = ^ — iu, f] = r] + i^y, and 

1 - (Mi/Ma)^ - 2i{Mi/M2) sin 47 



coth 2r] 



1 + (Mi/M2)2 - 2(Mi/M2) cos 47 
= ^R + t^i. (12) 



Note the non-trivial contribution from tan2ct; in Eq. (|l3). To diagonalize the symmet- 
ric and complex mass matrix, U~^m^U* , as is detailed in the next section, it is necessary 
to multiply the mass matrix on either side by the same phase matrix. This phase matrix 
is precisely e"*"^ 

Eq. (p!oD shows that, when m£, is hierarchical (^ ^ 1), the neutrino mixing angle 6 is 
small (tan^ ~ e~^^ ~ mi/ 1712) for generic values of the other parameters, /5, 77, and 7. 
However, when the denominator in Eq. ([10|) vanishes, 9 is maximal. This is the resonance 
behavior mentioned before. In general, the seesaw mechanism suppresses the neutrino 
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mixing angle. But when the resonance condition is met, it is enhanced and becomes 
maximal. 

This behavior is quantified in Fig.l, which is a 3D plot of the region sin^ 20 > 0.5, 
within the parameter space spanned by cos 2/5, 7 and M^jM^. This region consists roughly 
of two parts. One runs along the edge cos 2/3 ~ 1 and M1/M2 <^ 1, but 7 can take values 
between and 7r/4. The other region is tube-like, and "hugs" the back wall, 7 ^ 7r/4, with 
cos 2/? PS tanh2r7. It is striking how small the region for sir? 20 > 0.5 is. Outside of this 
region, which consists of most of the parameter space, sin^ 20 is tiny (~ (mi/m2)^). This 
result is the analog of the familiar focusing mechanism in relativity. When a relativistic 
particle decays, most of the decay products are contained in a forward cone of opening 
angle < I/70, where 70 = — This corresponds to the seesaw problem with 

the identification 70 = cosh 2^ ~ | (777-2 /mi). 

In Fig. 2, we blow up the region with a fixed cos 2/? 1. It is seen that there is 
considerable structure when sin^ 20 is maximal. In particular, the dependence on 7 is 
highly non-trivial. From the scale in the figure, we see that large values of sin^ 20 are 
confined in a very narrow region with width ~ [mxjm'2)'^ . Note also that, outside of 
the maximal sin^ 26' region, near the upper left edge of Fig.l (M1/M2 — > 0, cos 2/3 — > 1), 
sin^ 20 remains large. This region is characterized by extreme hierarchy of the Majorana 
masses {{rni/m2)'^ > (M1/M2) — > 0) and very small f5 ((1 — cos 2/3) ~ (7771/7712)^). 

Fig. 3 shows the contents of Fig.l in a 2D parameter space, with cos 2/3 = 1/2. It 
exhibits clearly the behavior of sin^ 20 near 7 = 7r/4. Here, the maximum of sin^ 20 is 
attained at cos 2/3 = tanh 2?] ■ tanh 2^ with 7 = 7r/4. Away from these values, sin^ 20 drops 
off quickly. The width of the peak is of order (7771/7712) in either A(Mi/M2) or A7. 

The behavior of physical neutrino mass ratio is depicted in Fig.4, which exhibits 
the region of near degeneracy, yUi//U2 > 0.5. We have chosen a log scale for M1/M2 to 
highlight the detailed structure near the upper left edge. A comparison with Figs. 1 and 

2 reveals the complementary nature of the regions of maximal sin^ 20 versus hierarchical 
/ii//i2. In the small region where sin^ 20 ~ 1, one also has fii/fi2 — 1- However, near the 
upper left edge, for (7771/7712)^ > M1/M2 — * 0, /ii//i2 can be small and at the same time 
sin^ 20 is large. For generic parameters, sin^ 20 ~ (7/21/7712)^, but the masses are also very 
hierarchical, {^1/^2) ~ (7771/7772)^. 

3 General properties of mass matrices 

In order to gain some insights about the results presented in the previous section, it is 
useful to study the general properties of symmetric, complex, matrices. We will first 
discuss the relations between different parametrizations of the mass matrices. These 
relations will shed light on the special properties of matrices of the seesaw type. They 
will also enable one to have a qualitative understanding of the results presented in Sec. ^ 
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3.1 Parametrization of neutrino mass matrices 

Within the framework of the seesaw model, the neutrinos are Majorana in nature, so that 
their matrices are symmetric and complex, in general. We first consider the case of two 
flavors, 

^-(bc)^ (13) 

Here, A, B, and C are arbitrary complex numbers. Without loss of generality, we assume 
that is normalized so that detN = +1, 

AC-B^ = 1. (14) 

This matrix can be diagonalized by a unitary matrix U, 

In terms of the eigenvalues (ni,n2), r] = | ln(n2/ni). A convenient choice for U is in the 
Euler parametrization 

The relation between the two parametrizations of is given by 



A B \ _ ( e2^'^(ch2r/ + C2/3 ■ sh2f]) S2/3 ■ sh2r/ 

B C I ~ \ 82/3- sh2r/ e-2*"(ch2r/ - C2/3 ■ sh2r/) 



(17) 



where we have used the notation t] = t] + i'y, ch2?7 = cosh 27], S2/3 = sin 2/3, etc. 

Note that because of the condition AC — B^ = 1, there are exactly four parameters 
in the three complex numbers A, B, and C. To understand the role played by the phase 
a, let us write 

N = ^{A + C) + i(A - C)as + Ba,. (18) 

The diagonalization of A^ is easy provided that the phase of A — C and B are the same. 
In general, we can multiply A^ on either side by the same phase, e"'""^^, 

^-ia^3j^^-iaas ^ l(e-2^"A + c^^^C) + ^(c'^^^A - e^'"C)a3 + Bai. (19) 

We now choose a so that the phase of {e~'^^°'A — e^*°C) coincides with that of B: 

arg B = arg(e-2'" A - e^'^C) . (20) 

In this case, the matrix e-*"o"3 ^Yg-^acs qq^^ be diagonalized by e'^^"'^ {e~^°"^^ Ne~^°"^^)e'''^"'^ , 
with 

2B 

tan2/5=^— -— = real. (21) 

Thus, given an arbitrary A^, we need first to determine the phase a by Eq. (^^. After 
which /3 is fixed by Eq. (|2T|), and then 7 can be read off from the diagonal matrix. 
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Note that from Eq.(|T]), 

Bie^'^'A* - e-^'"C*) = B*{e-'^'''A - e^'"C) . (22) 
Thus, Eg. (pop is equivalent to 

2a = aigiAB* + BC*) . (23) 
In addition, we may use Eq.(p^ in Eq.(pTD to obtain 

2\AB* + BC*\ 

tan2/3= • (24) 

Note also that, if \A\ = \C\, Eq.(pOD can not be used to solve for a, but Eqs.(p3D and ( ^4]) 
are still valid. The complex eigenvalues of Eq. (0) can be obtained from Eq. {\l7\} . They 
are given by 

6+2*? = Ae-^"" + B tan (3 (25) 
= -5tan/3 . (26) 

It is also useful to introduce another variable, 

C = IhA/C). (27) 

Using (, Eq. (|21| ) can be written as 

,,,2P^S^. (28) 



Also, from Eq. (0), 



sh(C - 2ia) 



ch2f] = ^(e-2^"A + e2^"C) 



Similarly, 
We thus have 



= Vac ■ch{C-2ia). (29) 

C2p-sh2f] = Vac ■sh{C-2ia). (30) 
C2/3 ■ tanh2r7 = tanh(^ — 2ia). (31) 



This relation can be regarded as a consistency check on the properties of A^. For instance, 
if C2/3 = (maximal mixing), it implies that Im^ = 2a, and that Re^ = 0. Another 
constraint is that the phase of tanh2?7 must be the same as that of tanh((^ — 2ia). 
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3.2 The seesaw transformation 




In the seesaw model, m,^ = m£)M^^mJ), it turns out that the properties of rrii, is closely 
related to M^^, when we choose a basis in which m^i is diagonal. We shall call the change 
from M^^ to m,^ a "seesaw transformation" (ST). In terms of the notation of the previous 
section, we define a ST from to a new matrix M by 

M = 

/ A' R' \ 

(32) 

It is seen immediately that B and AC are invariant {B' = B, A'C = AC), while 

A'/C = e-^^{A/C), (33) 

or 

C' = C-2e, (34) 
where (' = ^In^A'/C). If we assume that 

M = We^^^^W^, (35) 

we can use the results above to derive simple relations between the parameters pertaining 
to M and to A^. Thus, from the invariance of B under ST, we have immediately 

S2/3 ■ sh2f/ = S261 ■ sh2A, (37) 

i.e., S26i • sh2A is an invariant, independent of ^. One of its consequences is that the phase 
of A is tied to that of fj, since (3 and 6 are both real. In fact, if A = A + then 

tan 2(j) ■ coth 2A = tan 27 ■ coth 2// = constant, (38) 

independent of ^. In particular, if f/ = rj + in/A, sh2f/ is purely imaginary, then the 
imaginary part of A must also be 7r/4, i.e., the mass eigenvalues must have opposite signs. 
Moreover, given (3 and f], the relation exhibits the complementary nature of 6 and A, 
large 6 correlates with small A, and vice versa. This behavior was already discussed in 
connection with the results of Fig. 4 in Sec. ^ From Eq. (^), the invariance of Bj\fAC 
yields 

tan W = tan 2(3 , ^ ~ ^'"^ (39) 
^sh(C-2e-2zcu') ^ ^ 

When the ST is hierarchical, ^ ^ 1, it is clear that, for generic (, the angle 6 is suppressed 
(~ l/sh2,^ ~ (mi/m2)). However, if ^ ^ 2^, and if the phases in the denominator of 
Eq. ( |5P| ) cancel, then 6 becomes maximal. This was the behavior shown in Fig. 1. It 
should also be mentioned that Eq. (|39D reduces to Eq. ([lOD when one uses Eq. (|3TD. 
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As another application, we note that a quahtative understanding of Fig.l can be 
gleaned from Eqs. {^). Using Eqs. (0) and we have 

with u! = uj' — a given by Eq.(^). For ^ ^ 1, a necessary condition for large 6 is that 
|C| ~ (more precisely, \C\ < e-^«|v4| and \C\ < e'^^Bl). From Eq. (|13) , this means that 
ch.2fj ~ cos 2/5 sh277. However, since ch2f/ and sh2?7 have different phases, this equation 
has only special solutions. They are 1) 7 = 7r/4, so that both ch2f^ and sh2?7 are purely 
imaginary, and cos 2/5 ~ coth(2r7 + z7r/2) = tanh2r7. This last equation describes the 
shaded trajectory on the 7 = 7r/4 wall in Fig.[^. Another solution is 2) rj 00, so that 
ch2r/ ~ sh2r/ ~ e^'^e'^'>/2. Then ch2r/ - cos 2/? sh2r/ ~ e'^'^e^'^{l - cos2/5)/2. Since for 
?7 — >• 00, \B\ ~ sin 2/5 e^'', 9 can be large provided that e^^(l — cos 2/5) -C sin 2/5. This 
solution corresponds to the shaded region in Fig.[I| with 7 7^ 7r/4. 

In summary, the neutrino mixing angle 6 can only be large if the (2, 2) element of 
M^^ is small, |C| ~ 0. The precise value depends on phases and possible cancellation 
between A' and C. Note that, in the literature, a number of studies has concentrated 
on the case of M^^ being a real matrix. For large mixing, two types of M^^ have been 

identified. 1) M^^ ~ {l e^^ )' " ( 1 ) ' °^ ~ ^ ( 1 ) ' 

detM^^ = +1. These are special cases of |C| ^ 0, corresponding to the two end points of 
the shaded region in Fig.!, with coordinates (Mi /M2, cos 2/5, 7) ~ (0,1,0) and (l,0,7r/4), 
respectively. Our analysis shows that care must be taken when we have small deviations 
from these forms, which arise naturally in models constructed from a presumed broken 
symmetry. The narrowness of the shaded region means that a viable solution can be easily 
thrown off course by small perturbations. An example of such sensitivities is known in 
the renormalization group running effects, which we will discuss in the next section. 



3.3 Renormalization 

It turns out that our general analysis has an immediate application to the renormalization 
group equation (RGE) analysis of the neutrino mass matrix. We briefly comment on this 
connection. A full account will be given elsewhere^. 

In the SM and MSSM, the RGE running of the neutrino mass matrix has been very 
extensively studied[§, ^. For simplicity, we only consider the two flavor problem with 
(i/^, Ut). The RGE for the effective neutrino mass matrix is given by, 

^rrii, = -(nm^ + m,jP + P'^m^) , (41) 
at 

where k is related to coupling constants, t = In fi/Mx, and to a good approximation, 

P ~ ^ X(l - ^3) , (42) 
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where x is given by in the SM and —y^jl in the MSSM, with y^. and being the r 
Yukawa couphng in the SM and MSSM respectively. The solution to RGE is 

m^it) = e~'''*e«"^m^(0)e«"^ , (43) 

where k' = k + 2x, ^ = X^i and we have ignored the t-dependence of the coupling constant 
so that / ndt ~ nt, etc. 

It is convenient to factor out the determinant 

rriy = ^mim2M . (44) 

Then , 

^Jmi{t)m2{t) = e-"'ymi(0)m2(0) (45) 
M{t) = e«'^W(0)e«'"^ (46) 

Thus, while the overall scale ^m\m2 has a simple exponential dependence on t, the 
running of M, which contains the mass ratio and the mixing angle, is just a seesaw 
transformation defined in the previous section. The difference from the traditional seesaw 
model is that, instead of ^ ^ 1, for the RGE running ^ is usually small (~ 10^^ in the 
SM). Nevertheless, the exact and analytic formulae given in Eqs.(^|l^) are valid solutions 
of the RGE. Detailed analysis of their properties will be given in a separate paper. We 
only note that, according to Eq.(^), there is a (complex) RGE invariant, 

sin 2e{t) sinh 2A(t) = sin 2^(0) sinh 2A(0) , (47) 

where Q and e^^ (A = A + z0) are, respectively, the physical neutrino mixing angle and the 
mass ratio. This equation can be used to determine A(t) and 0(t), once Qii) is obtained 
from Eq.(|lOD. We should also emphasize that, because the solution can exhibit resonant 
behavior, large effect can result even for very small running (^ -C 1). 



3.4 Parametrization of the three flavor matrix 

As is clear from the previous discussions, the Euler parametrization is the most convenient 
for dealing with the two flavor problem. The generalization to three flavor, then, amounts 
to parametrizing an ^[/(S) element in the form (phase) (rotation) (phase). However, there 
are altogether eight parameters in ^[/(S) while each phase matrix can only accommodate 
two. So there must also be an additional phase matrix contained in the rotational part 
of a general S'[/(3) matrix. This decomposition is of course none other than the familiar 
CKM matrix decomposition. Thus, for three flavors, the analog of Eq. (^) is 

_ g«{e3A3+£8A8)gie7A7g«e5A5g«<53A3gie2A2gi(£3-'*3+e8-'^8) ^^g'j 

Like the CKM representation, the phase factor e**^^^^ could be put in a different location, 
or one could use another diagonal A matrix. 
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The seesaw problem for three flavors again aims at rewriting the matrices so that m,^ 
is given as in Eq. (^), with W assuming the form of Eq. (^). As in the two flavor case, the 
exterior phase factors of W do not contribute to neutrino oscillations. An exact solution 
for the three flavor problem, however, is not easily obtained owing to the complexity 
of computing finite matrices involving the A matrices. In Sec. ^, we will present an 
approximate solution to the three neutrino problem. 



4 A unified approach to fermion mass matrices 

Our general analysis of the properties of the seesaw model suggests a unified picture of the 
quark and neutrino mass matrices. As was discussed in Sec. 0, the physical mixing angle 
of a seesaw model is quite small, in general, but can be maximal when special conditions 
are met. We will now present arguments which can associate these regions to the quark 
and neutrino mass matrices, respectively. For simplicity, our discussions are restricted to 
the case of two flavor only. 



4.1 Quark mass matrices 

It has been known for a long time that quark mass matrices can be adequately described 



in a seesaw form|10, 11 



m 



Hi 




(49) 



where a,b,c are arbitrary complex numbers, all of the same order, and H2/f^i ^ 1- Here 
we have used the arbitrariness in m to demand that it be complex and symmetric, in 
contrast to the usual choice that m is hermitian. 



If the physical masses are denoted as mi and m2, then Eq. (|49D implies that, for generic 
values of a, b, c, 

mi/m2 ~ (/ii//i2)^ • (50) 



While the mixing angle satisfies the well-known relation|10 



sin^ 9 ~ nil/ 1712 ■ 



(51) 



This result was derived first for real matrices, and remains valid for complex case, as 
discussed in Sec. H. It has served as a model for quark mass matrices for a long time. 

Physically, a symmetric and complex mass matrix can be derived by a symmetry 
argument. Since the mass term in the lagrangian is given by g^Mg^, a symmetric mass 
matrix can be naturally obtained by imposing a discrete Z2 symmetry: 



Ql ^ Qr 



(52) 
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If we further impose a gauged horizontal symmetry, such as a f/ (1) symmetry a la Froggatt 



and Nielsen [|r^, then we are led to a mass matrix in the form of Eg. . For instance, 
we may take the horizontal charge assignments (0, 1) for (g^i, 5^2) and (1, 0) for {qm, gm)- 
The charge assignments for the mass matrix qiMgn is 

Qm-(^1 ly (53) 
The Froggatt-Nielsen mechanism then calls for a mass matrix of a form 

M^^t), (54) 



as in Eq. 

4.2 Neutrino mass matrices 



In Sec. we have found that, in order to have a large physical mixing angle, it is 
necessary that the (2,2) element of be small, i.e., \C\ ^ for M]^^ ~ ^ 5 C 

Since we have used the normalization detM^^ = +1, the Majorana mass matrix is given 
by 

= ( -B :f ) ■ (55) 

The condition |C| ~ simply means that is itself of the seesaw form. The condition 
\C\ ~ is not sufficient, however, to guarantee a large mixing angle, which is a consequence 
of further constraints on Mj^. We will not attempt a detailed model construction here. 
We only note that, as emphasized in Sees. ^ and ^, the mixing angle is very sensitive to 
small variation of the parameters in M^. In particular, if a model is based on symmetry 
arguments, symmetry breaking effects have to be weighed carefully. 

In summary, both the quark and neutrino mass matrices can be adequately described 
in the seesaw form. Their difference arises from the Majorana sector, which is itself of the 
seesaw form. This last requirement can lead to large mixing in the effective neutrino mass 
matrix. The sensitivity to small changes in the parameters calls for a careful examination 
which should also include three flavor effects. Detailed model construction along these 
lines will be attempted in the future. 



An approximate solution to the three flavor prob- 
lem 



In Sec. |3]^, it was pointed out that the three flavor seesaw |T2[ problem amounts to 
rearranging products of matrices in SL{3,C). Since a general, analytical, solution is not 
available, we will turn to an approximate solution which is physically relevant. 
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For the three neutrino problem, it is known that the (23) angle is near maximal, the 
(13) angle is small, and that the (12) angle is probably large. This suggests that, to a 
good approximation, the three flavor problem can be decomposed into two, two flavor 
problem. To implement this scenario, let us consider the 3x3 matrix M^^, 



Mr' 



j A B D\ 

B C E 
\D E F J 



(56) 



The neutrino matrix, with m^) diagonal and U = I for simplicity of presentation, since 
the general case can be easily incorporated as in Eq. (0), is given by 



m,. 



( ^1 



\ 



1712 




m2 



(57) 



ma 



It is convenient to introduce, in addition to the Gell-Mann A matrices, Ag and Aio, 





(5^ 



\/3A 



10 



(59) 



We may now write 

/ nil 

V 



m2 



nil 



rris 



C = o In(m3/m2). 



Then, the (23) submatrix of m^, can be diagonalized. 





M 


B 


^\ ( 




B 


C 


E 1 e"^^'' = U 1 




{ D 


E 





A B' D' 



B' A' 
D' 







u 







( A 


B 




B' 


dI 




{ D' 



-5A9 



(60) 
(61) 

(62) 

(63) 
(64) 



and (A', S') are the eigenvalues. Although we could have chosen a proper normalizing 
factor so that the (23) submatrix has det = ±1, as in Sec. 0, for this problem it is simpler 
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not to do this and A'S' 7^ 1, in general. If we absorb e*'^'^^ by defining the new variables 
(S",L'",A",E")= (e*'^S',e-*Ti:'',e2*'^A',e-2'TE'), and since A7 and Ag commute with the 
remaining Dirac matrix, Eq. (|57D becomes 



X 



A 


B" 


D" 


B" 


A" 





D" 







mi 







X 



T 



(65) 



(66) 



Now, the (13) rotation is controlled by However, we must first make sure 

that they have the same phase (with the approximation |mfy4| <C 1 777.2/713 S"|). To this end 
let us multiply 777.,^ by e*'^^'^"' on either side, and choose uo so that e^^D" and e^*'^S" have 
the same phase, arg(e~*'^D") = arg(e^*'^S"). In this case, we can rotate away the (13) 
element of 777,^ without changing its other elements by assuming that the angle of rotation 
is small, |777iD"| <C y^m^m^S'' . We have (with A" = A — D"'^ approximately. 





/ N' 


B" 







B" 


A" 







[ 








(67) 



Y = g*"-^9g-*^(v^-^lo)g«/3-^7g-«'/'A5gia;(v^Aio) 



777i 



tan 7/^ = (777ie-*'^D")/(y^«i^e'"^S") = real. 



(68) 
(69) 



After this somewhat laborious route, we see that the diagonalization of my can be 
finally achieved by working solely in the (12) sector. The crucial assumption for the 
success of this procedure is that ianip <^ 1. Otherwise the (13) rotation e'^"^^'^ will 
generate non- negligible elements all over the matrix m^. Although the exact condition 
for tan 7/; <C 1 seems complicated, in practice, as long as the elements B and D in 
are reasonably small, the approximation is valid. 

Fortunately, it is known that in reality the physical (13) rotation angle is small. This 
means that for any successful m^, the above approximation is appropriate. In this case, 
the three neutrino problem is reduced to two, two-flavor problem. In particular, two 
popular scenarios, the bimaximal or single maximal models, can be accommodated. 



6 Conclusion 



Recent experimental data have revealed two striking features of the intrinsic properties 
of the neutrinos. One, as expected, they are very light. Two, perhaps surprisingly, at 
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least some of their mixing angles are large, or even maximal. The seesaw model provides 
a natural explanation of the lightness. However, the story of the mixing angles is more 
complicated. In the seesaw model, the neutrino mixing matrix can be written as UW, 
where U comes from the left-handed rotation which diagonalizes the Dirac mass matrix, 
and W, defined in Eq. (|^), is induced from the right-handed sector of the model. For two 
flavors, the analytic solution for W shows that, when there is a mass hierarchy in mu, 
the mixing angle in W is greatly suppressed for most of the available parameter space. 
However, in a very small region, which we exhibited explicitly in Sec. |^, the mixing angle 
can be large. In addition, this region may be divided roughly into two parts. In one, 
characterized by 7 ~ 71/4, the physical neutrino masses are nearly degenerate. In the 
other, in which the Majorana mass eigenvalues are hierarchical, the neutrino masses can 
be either hierarchical or nearly degenerate. This behavior of W has interesting theoretical 
implications. 

Since the neutrino mixing matrix is given by UW, there are three obvious possibilities 
which can lead to large mixing. A) U contains large angles but ~ /; B) both U 
and W contribute appreciably and they add up to form large mixing; C) f/ ~ / but the 
large angle is in W. Corresponding to these possibilities we have three different physical 
scenarios. A) With 1^ ~ J, the physical neutrino masses are highly hierarchical. The 
burden for the model builders is to find a credible theory which makes U almost maximal 
naturally. B) This scenario seems the least likely to be implemented. This is a "just-so" 
solution whereby the Dirac and Majorana sectors must conspire to make the resultant 
angle large. C) Here, ?7 ~ J is quite reasonable from quark-lepton symmetry, which leads 
naturally to f/ ~ Uckm- The challenge is to find a mechanism whereby the parameters 
in the seesaw model lies naturally in the narrow range for large mixing. 



In Sec. [4.2| , we have identified a necessary condition for large mixing, namely, that 
the Majorana mass matrix is also of the seesaw type. This result suggests a universal 
seesaw mechanism for both the quark and neutrino mass matrices. The quarks can take 
advantage of the general solution, resulting in small mixing and hierarchical masses. For 
neutrinos, the existence of Mr can then lead to large mixing angles. More detailed studies 
are necessary to implement this scenario. 

Our general analysis of symmetric and complex matrices also has an immediate appli- 
cation to the RGE running of the neutrino mass matrices. Exact and analytic solutions 
of the RGE are found, in addition to a (complex) RGE invariant which relates explicitly 
the running of the mixing angle, the mass ratio and its phase. 

The analyses given above are for the case of two flavors. However, in the approximation 
of a small (13) angle, we have found that the three flavor problem is reduced to two, two 
flavor problems. We thus do not expect qualitatively different physics for this case. 

In conclusion, the neutrino mixing matrix (masses and mixing angles) implied by the 
seesaw model has rather intriguing properties. To accommodate large mixing angles, 
there are just a few limited options available. These conditions should be helpful in the 
search of a viable neutrino mass matrix. We hope to return to this topic in the future. 
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Figure 1: Region in which sin^ 261 > 0.5, with tanh2^ = 0.9998, or (mi/ms) = 0.01. Note 
the log scale used for (M1/M2). 
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Figure 2: A plot of sin^ 2^ vs. (M1/M2) and 7, with cos 2/3 = 0.9999, tanh2^ = 0.9998. 
Note the expanded scale of (M1/M2). 
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Figure 3: Typical behavior of sin^ 26' for cos 2/5 < tanh2^. Here, cos 2/3 = 0.5, tanh2^ = 
0.9998. 
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Figure 4: Region in which the physical neutrino masses are nearly degenerate, with 
A*i/a«2 > 0.5, tanh2^ = 0.9998. 
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